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'st ■ Abstract 

The competing influence of two types of interactions on surface diffusion is investigated: one of 

Q , them, e, acts between adsorbed particles, while the other one, r/, between activated and adsorbed 

^ ' particles. To this end, a specific lattice-gas model on a triangular, square, and hexagonal lattice is 

"t^ ' considered with an attractive e and an attractive or repulsive r], both restricted to nearest neighbors. 

-)— » \ 

"^ , For all three lattices the infiuence is qualitatively the same. Namely, when 77 is neglected, then e 

-)— > 

rH \ is shown to decelerate diffusion with a rate exponential in e and inverse temperature. Moreover, 

F^ ' when r/ is present and e is fixed, then a sufficiently attractive r/ (relative to e) accelerates diffusion 

Q . exponentially fast in r/ and inverse temperature. However, quite surprisingly, a repulsive or slightly 

o ■ 

attractive r] has practically no effect on diffusion. Finally, when r] is set proportional to e via 
Kj^ a parameter a, surface diffusion is exponentially accelerated (decelerated) for a above (below) a 

\^ • threshold value equal to 4, 5/4, and 7/8 for the three lattices. Thus, an r/ of strength comparable 

^: 

\^ , to e is enough to boost diffusion on the square and hexagonal lattices, while an r/ of a rather large 



o 



X 



strength is needed for the boost on the triangular lattice. 
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I. INTRODUCTION 

The importance of surface diffusion as a mechanism of mass transport has been ac- 
knowledged in various fields, such as chromatography,- sintering,-*^ microelectronics,-!^ civil 
engineering,-'^ heterogeneous catalysis,- or neuroscience.-ii^ A widely used approach to sim- 
ulate surface diffusion at a microscopic level is to employ lattice-gas models.— "^^ In these 
models the diffusion is given by a potential relief of the surface. While most of the time the 
particles stay at adsorption sites (minima of the relief), occasionally they perform random 
jumps to adjacent vacant sites. If the jumps are fast enough, the microstates of the system 
can be represented by occupation numbers (one number for each site), as in a lattice gas. 

Interactions between particles may strongly influence surface diffusion. Indeed, attractive 
interactions between particles adsorbed on the surface restrain their migration and, thus, 
slow down diffusion. By the same token, repulsion between adsorbed particles boosts diffu- 
sion. This basic picture is obscured by the fact that particles in the activated state — those 
performing jumps between two sites and located at saddle points of the potential relief — 
also interact with the adsorbed particles within their vicinity, which may have a significant 
impact on diffusion as well. However, the effect of this interaction is opposite: if it is at- 
tractive, it tends to lower the effective barriers of the jumps, which leads to higher jump 
probabilities and a boost in diffusion; if it is repulsive, it should slow diffusion down. Us- 
ing Monte Carlo simulations at temperatures above the critical point, the effect of each of 
these two types of interactions on diffusion was already studied separately, corroborating 
the above-described anticipated behavior.— i^^ A combined effect of both types of interaction 
has not been studied, though. Hence, it still remains unclear which effect is dominant and 
under what conditions. This point is investigated in detail in the present paper. 

To this end, we employ a simple lattice-gas model with an attractive nearest-neighbor in- 
teraction, e, between adsorbed particles, and a different (attractive or repulsive) interaction, 
rj, between an activated particle and its closest adsorbed particles. A triangular, square, 
and hexagonal lattice will be considered, all of which will turn out to exhibit the same type 
of behavior in surface diffusion coefficients. 

Our analysis is based on two basic assumptions. First, we will assume that the surface 
coverage varies only very slowly with time and space, i.e., that the local equilibrium ap- 
proximation is applicable. Then the chemical and jump diffusion coefficients, Dc and Dj, 



can be expressed via purely thermodynamic quantities. As various comparisons with kinetic 
simulations have revealed, the results obtained within this approximation are reliable even 
farther away from equilibrium,—"— indicating that the approximation is actually quite ro- 
bust. Second, we assume that the temperature is low (sufficiently below the critical point) 
so that ordered phases of adsorbed particles occur in the system. Then explicit formulas for 
coverage dependences of D^ and Dj are available,— and the needed thermodynamic quan- 
tities can be evaluated (to any desired precision, in principle) via cluster expansions. In 
particular, one can evaluate various multi-site correlation functions that inevitably arise due 
to the presence of an activated-state interaction,^- Previously, the correlations were obtained 
only by approximative techniques, like the cluster variation method^is^S or approximations 
via two- or three-site correlations.— ^^^ Rather crudely, the activated-state interactions were 
often neglected to work only with the correlations between nearest-neighbor sites.— i^^i^^i"— 
The paper structure is as follows. In Section [TTl we introduce the studied lattice-gas 
model and give expressions for the coefficients Df. and Dj in terms of the surface coverage, 
isothermal susceptibility, and a correlation factor. In Section IIIII we then investigate the 
effects of the two interactions e and r] on these coefficients, employing general formulas for 
the coverage dependence of Df. and Dj. Upon evaluating the quantities involved in these 
formulas for our model, we obtain explicit dependences of Dc and Dj on the interactions, 
allowing us to analyze in detail their influence on surface diffusion. A summary of our results 
and concluding remarks are given in a final section. 

II. MODEL OF SURFACE DIFFUSION 

We shall consider an ideal solid surface whose potential relief minima form a triangular, 
square, or hexagonal lattice. Each lattice site, x, is either vacant (an occupation number 
Vx = 0) or occupied by an adparticle {u^ = 1). An interaction between adparticles is assumed 
to occur only between nearest neighbors (nns), the corresponding interaction energy being 
attractive, e < 0. The model Hamiltonian is given as 

H = eN2-fim, (1) 

where Ni {N2) is the number of occupied sites (occupied nn pairs) in a given microstate and 
fi is the chemical potential. 
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FIG. 1: The three lattices with a jump between a pair of nn sites 1 and 2. The segments 5" are 
indicated by shaded areas. 

The model is equivalent to the Ising model of ferromagnetism and was already used to 
simulate surface diffusion at high temperatures T.— "— >^ In particular, it was employed to 
study the effects of an attractive interaction e on surface diffusion on a triangular lattice (for 
T > 0.94|£:|/A;^) and on a square lattice (for T > 0.568\6\/kB and e attractive).— 'i^ Note 
that these temperatures are above the critical values 0.9102\e\/kB and 0.5673\e\/kB for the 
two lattices (for a hexagonal lattice the critical point is 0.3797|£:|/A;b)-— 

At low temperatures model ([I]) is known to have only two phases: a fully vacant phase 
at fi < fit and a fully occupied phase at fi > fit; at fi = fit both phases coexist.—"^ Thus, 
the models exhibits a first-order transition at fit between the fully vacant and fully occupied 
phase. The point fit = qe/^, where q is the lattice coordination number (equal to 6, 4, and 
3 for the three lattices). 

Migration due to surface diffusion is represented by jumps of adparticles to nearby vacant 
sites. For simplicity, we shall assume that these jumps are uncorrelated and occur only 
between nn sites. The influence of adparticles on the activation energy of a jump will be 
taken into account via an (attractive or repulsive) interaction, 77, that is different from e and 
acts between an activated particle and its closest adparticles.—*^ Thus, this interaction is 
associated with segments, S, of more than two adsorption sites (see Fig. [1]): a pair of nn 
sites between which a particle jump is performed, plus the closest sites to the pair's center 
(over which a saddle point of the potential relief is located). 

Relevant many-particle transport parameters for surface diffusion are the chemical and 
jump diffusion coefficients, D^ and Dj, respectively. The former coefficient describes surface 
mass-transport processes of the system and is defined via the Fick's first law; the latter 
describes the asymptotic behavior of the mean square displacement of the system center 
of mass. They are related by the Kubo-Green equation, Dc = {/36/x)Dj,—'^ where 6 



is the surface coverage, x is the isothermal susceptibihty, and /3 = l/ksT is the inverse 
temperature. 

As long as the coverage varies only very slowly with time and space (i.e., the local equi- 
librium limit is applicable), purely thermodynamic quantities are sufficient to approximate 
the two diffusion coefficients,— '^^'2^'^"— 

D, ^ I^oe^^-^, Dj ^ Doe^"^, (2) 

where Dq is the diffusion coefficient of non-interacting particles. 

The correlation factor P is associated with the interaction of activated particles and may 
be written asM-^ 

P = ^e-("s-"*^)''Vf. (3) 

M 

The summation is over all subsets M of a segment 5* that contain the pair of nn sites between 
which a jump is performed, including S itself (see Fig. [2]), and um > 2 is the number of 
sites in M. Moreover, pm is the statistical average that the sites of M are vacant and the 
remaining sites in S are occupied. For example, for a triangular lattice and M = M3 we 
have ns = 4, um^ = 3, and pms = ((1 — ^i)(l — '^2)(1 — t'3)'^4), where 1, ... ,4 are the sites 
of S and 1 and 2 are the nn sites associated with a particle jump. Taking into account 
that there may be several sets M yielding the same value of the average pm, one explicitly 
hasi^^ 

P = e-^^PpM, + 2e~^PpM, + PS (4a) 

for the triangular lattice and 

P = e-^''V/.+4e-3'"3pM3 

+ 2e-=^''V + 4e-''Vf5+P5 (4b) 

for the square and hexagonal lattices, where the shorthand p4 = pu^^ + PM42 + VMts ■ 

III. RESULTS AND DISCUSSION 

According to Eq. ([2]), in the local equilibrium limit the diffusion coefficients D^ and Dj 
become thermodynamic quantities. A rather powerful statistical mechanical machinery2i~i22 
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FIG. 2: The statistical averages pM of finding holes at the sites of M (indicated by circles) and 
particles at the other sites of S (indicated by disks) for the three lattices. A given pM may 
correspond to several sets M of which only one is shown (for example, there are two sets M3 on a 
triangular lattice associated with pm^)- 

was developed to study such quantities at low temperatures. We shall apply it here to 
investigate the coverage dependences of Dc and Dj, focusing on their variations with the 
interactions e and 1]. Before doing so, let us mention two peculiar features of systems at low 
temperatures that underlie the forthcoming results. 

The first one is the structure of phases. At low temperatures a typical microstate in a 
given phase looks as a 'sea' of a ground state in which isolated 'islands' of non-ground-state 
configurations are scattered.—"— For the fully vacant (occupied) phase in model ([1]) this 
corresponds to a sea of holes (particles) with islands of particles (holes). The islands are 
small, with diameters of order not exceeding InA^, yet their bulk density is non-zero, of 
order exp(— c/3) <^ 1. In fact, the smallest islands (of one or few sites) yield the dominant 
contributions to the system's free energy. As a result, in a given phase the coverage has 
a practically constant value corresponding to the ground state associated with the phase 
(a value or 1 for model ([1])), the deviations being just of order exp(— c/3). Hence, at low 
temperatures a significant change in the coverage can possibly occur only at or very near a 
transition between phases. 

The second feature is the smoothing of a discontinuity in a thermodynamic quantity (such 
as the coverage) at a first-order phase transition in a finite system. At low temperatures the 
profile of the smoothing is identical (given by the function tanh) for a large group of lattice- 
gas models.—'^ This allows one to described the dependence of thermodynamic quantities 



on the coverage in a unified way. 

A. Formulas for diffusion coefficients 

These facts can be used to obtain general low-temperature formulas for the coverage de- 
pendences of Dc and Dj associated with a first-order phase transition between two phases.— 
In the case of model ([T]) they read 

^ Dpe^^^ / Po P. \ 

"^ AON \e-eo e,-eJ' .^s 

where A^ is the total number of adsorption sites in the system, and the constants 9o, 0, and 
Po, P, are the coverages and correlation factors, respectively, evaluated at the transition point 
fit in the fully vacant (symbol 'o') and fully occupied (symbol '•') phase. The shorthand 
A^ = 6, ~ 60 is the coverage discontinuity at the transition. 

Formulas ([5]) work for a system with periodic boundary conditions and for coverages 
60 < < 0,- This is almost the whole interval 0<^<1, as^o~0 and ^, ^ 1 at low 
temperatures (see Eq. ([6]) below). The formulas show that the coefficient Dc behaves as 
a sum of two hyperbolas diverging at 0o and 0,, respectively, and the coefficient Dj as a 
hyperbola B + C/0 diverging at (see Fig [3]). No divergences in D^ or Dj actually occur, 
however, because for < ^ < ^o + 5 and 0, -\- 5 < < I with 5 ~ N^'^^^ <^ 1 one of 
the two phases prevails in the system, and formulas different from Eq. ([5]) become valid.— 
Nevertheless, in the following we shall use only Eq. ([5]) since it is applicable to practically 
all coverages. 

The interactions e and rj affect the diffusion coefficients via the parameters Po,P,, 0o,0», 
A0, and fit that we shall now evaluate. Note that rj appears only in the correlation factors 
Po,P, (via the exponential exp(— r^/?)), whereas the interaction e affects Po,P, (via the 
averages P°m,P*m) ^s well as the coverages 0o, 0, and location of the transition point fif 

B. Single-phase coverages and correlation factors 

At low temperatures, the leading contributions to the free energy in a given phase come 
from single-site thermal perturbations of the ground state. In the vacant phase such a 
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FIG. 3: The coverage dependence of the diffusion coefficients as given by Eq. dH) for various 
possibihties of parameters values. For the jump coefficient Dj ~ i? + C/6 the parameters B = 
P, — Po and C = PqO, — P,9o cannot be simultaneously negative. 

perturbation corresponds to microstate with an energy excess /S.H = — yU > over the vacant 
ground state, while in the occupied phase the excess is AH = fi — qe > 0. These excesses 
correspond to the Boltzmann factors exp(/i/3) and exp[(g£: — fi)f3], respectively. Hence,— i^ 
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(6) 



Note that this yields AG ^ 1 - 2exp{qe/3/2). 

In an analogous way, we may obtain the single-phase averages p^^ and p\.j at the transition. 
Consider, for example, the average PM42 '^^ ^ square lattice. It corresponds to addition of 
two nn adparticles to the fully vacant ground state (see Fig. |2]). Similarly, p'^^^ corresponds 
to removal of four adparticles in an elementary square from the fully occupied ground state. 
The energy excesses of these microstates are AH = e — 2^ and AH = 4yU — 125, respectively. 
Thus, within the leading-order approximation, p°m^^ ~ exp[(2yUt — £:)/3] = exp(3e/3) and 
P'ma2 ~ exp[(12£: — 4yUi)/3] = exp(4£:/3). In this way the leading terms in all averages p\j and 
p\,j can be deduced (see Table [T]). Although this evaluation is heuristic, it can be shown to be 
actually correct. ^° It must be taken into account, however, that there may be two or three 
different types of perturbations of a ground state corresponding to the leading term in p\j 
or p^, yielding a non-unit multiplicative prefactor in such a case. 

The single-phase correlation factors Po and P, are given by Eq. ([3]) with the averages pm 
replaced by p^j and p\j, respectively. Combined with Table [U explicit expressions for Pq 
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TABLE I: The leading-order terms in the sing 



transition point fi = Ht (adapted from Ref. 
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e-phase statistical averages p^j , p'j^j evaluated at the 



and P, readily follow. Namely, 



p^^g(6e-2,)/3^2e(3.-,,)/3^^^ 



for the triangular lattice, 



P. ?« e^^^"''''^'^ + 4e(^^-3'?)/^ + 4g{5e-2,y)/3 



for the square lattice, and 



Po ^ 4e(6e~4'))/3 _^ gg(9e/2-3»;)/3 

p. ^ e^^^"^'')^ + 4e(5^/2-3^)/3 + 3g(3e-2r,)/3 
^8g(7e/2-r,)/3^4g4./3 



(7a) 
(7b) 



(7c) 
(7d) 



(7e) 
(7f) 



for the hexagonal lattice. 
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FIG. 4: The dependence of the correlation factor (a) in the vacant phase and (b) in the occupied 
phase on the activated-state interaction r] for the three lattices (denoted as Tri, Sq, and Hex) and 
a fixed interaction e = — 3//3. The ranges in which the terms corresponding to the sets M = M2 
and M = S are dominant are indicated, and the sum of these two terms for each factor and each 
lattice is shown by a dotted line. 

The rj dependence of the correlation factors for a fixed e as given by Eq. ([7]) is shown in 
Fig. m Obviously, there is a single dominant term in Pq and P, for rj either well below or 
well above a certain value, rjo and rj,, respectively. Below these vales, the dominant term 
corresponds to the set M = M2, while above them it corresponds to the set M = S*. On the 
other hand, near the values rjo and rj, three or more (or even all) terms become essential. 
Note that t]o = 3e and rj, = —e for a triangular lattice, while to rjo = 'iej'l and r^, = — e/2 
for square and hexagonal lattices. Finally, comparing the factors Pq and P, at a given r^, we 
observe that Pq prevails at high r/, whereas P, prevails at low r^. 

C. Influence of interactions on diffusion coefficients 

Equations (E]) ^ (EI) provide explicit dependences of the diffusion coefficients D^ and Dj on 
the interactions e and r^. We may simplify these dependences and make thus their further 
analysis easier, if we use that for any t] only one term in Po and P, is dominant in the 
coefficients. 

To see this, we rewrite the two expressions in the parentheses in Eq. ^ as V /{9 — 9o){9, — 
6) and V/O, respectively, where V = {P.- Po){9-9*) + {P, + Po)A9/2 with 9* = {9o + 9,)/2. 
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From Eq. ([3]) and Table [T] we get 

M=M2,S 

^e--^^pl,^±p^, (8) 

where w = ns — nM2 is equal to the number of adparticles interacting with an activated 
particle (2 for the triangular lattice and 4 for the square and hexagonal lattices). Hence, only 
the terms corresponding to M = M2 and M = S are important in P, ± Po, while the other 
terms are exponentially suppressed (their relative contributions are of order exp{eP/2) ^ 1 
or less). Substituting Eq. ([8]) into V, we readily get that the parentheses in Eq. ([5]) can be 
approximated by the same expressions with Po and P. replaced by pg and exp{—wrj/3)pl^^, 
respectively. We thus conclude that the term p°g from Pq and the term exp{—wri/3)p'^^ from 
P, prevail in the diffusion coefficients at low temperatures. 

Combining this with p°^ ^ 1, p*^^ ^ exp[(g — l)£:/3], 60 ~ 0, and 6, ^ 1, we get 






Dj ^ Doe^'^/^^e^^"-'^'-^^^" + ^— ), 



with a possible exception of coverages close to and 1. Only two lattice parameters appear 
in these two formulas: q and w, i.e., the number of adparticles interacting with a particle in 
the adsorbed and activated state, respectively. Note that the term exp{—wri(3)pl^^ prevails in 
Dc and Dj when 77 < a*e with a* = {q— l)/w > 0, i.e., when the activated-state interaction 
is sufficiently attractive (more than 5/2, 3/4, and 1/2 of e for the three lattices). On the 
other hand, if 77 is repulsive or less attractive than a*e, then the term pg prevails. 

Equation iQ will allow us to easily analyze the influence of the interactions e and t] on 
the diffusion coefficients. We shall discuss three specific cases. 

Case 1: The interaction in the activated state is neglected [r] = 0) and only the influence 
of 6 occurs. Then Eq. ([9]) reduces to 

D.^^l Z,.«Z,„e-/^i^. (10) 

Thus, the diffusion coefficients decrease as e is more and more attractive (as is expected). 
The decrease is exponential and occurs primarily due to the chemical potential term 
exp(yUi/3) = exp{qe(3/2). This behavior of the diffusion coefficients is illustrated in Fig. |5l 
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FIG. 5: The coverage dependences of the diffusion coefficients for the three lattices when the 
activated-state interaction is neglected (77 = 0). The solid lines correspond to Eq. (j5]) and the 
dashed lines to the approximation from Eq. ([TO]) . The curves are labeled by the value of e/3. 

Note that the approximation from Eq. f lTOj) is in very good agreement with Eq. (jS]), faiUng 
only for Dc near 6 = 1. 

Remark. When the interaction rj is neglected, only the single-site and nn-site statistical 
averages occur in the correlation factor p2kd^^-2^ Indeed, for 77 = Eq. (jl]) reduces to 
P = 1 — 26+ {v\V2) with 6 = {vi). Then the free energy of model ([1]) is sufficient to evaluate 
P. However, as soon as rj is taken into account, multi-site averages pm inevitably arise, and 
the model's free energy cannot yield P. Instead, one should use the free energy of a more 
complex version of model ([T]) in which the needed multi-site interactions are added.— 

Case 2: The interaction e is fixed and the value of r] varies. Then it would be antic- 
ipated that the diffusion coefficients decrease with rj repulsive and grow with rj attractive. 
Nevertheless, this is true only when the term exp(— wr]/?)^*^^ is dominant, i.e., when rj suf- 
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FIG. 6: The coverage dependences of the diffusion coefficients for the three lattices when the 
adsorbed-state interaction is fixed (e = — 5//3) and the activated-state interaction ij varies. The 
solid lines correspond to Eq. ^ and the dashed lines to the approximation from Eq. Q. The 
curves are labeled by the value of ij. Above certain values of r/ the curves are practically identical. 

ficiently attractive, r] < a*e. Then the diffusion coefficients grow with rj exponentially fast 
as exp(— w?7/3). However, for r] repulsive or slightly attractive, t] > a*e, the term pg^l 
prevails, and the diffusion coefficients are almost rj independent, Dc ~ Doexp{fitP)/N6 and 
Dj ^ Dq exp{^tf3) {1 — 9) / 9 . Thus, no influence on surface diffusion occurs due to t], contrary 
to expectations, as is shown in Fig. O Again, the approximation of Eq. i^ is in very good 
agreement with Eq. (jS]), except near 9 = 0. 

Case 3. Finally, assume that the two interactions are mutually related. For simplicity, 
let them be proportional to each other, rj = ae, where a is positive or negative or zero. Then 
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FIG. 7: The coverage dependences of the diffusion coefficients for the three lattices when rj = ae 
with a = a* < ao- The sohd hnes correspond to Eq. ([5]) and the dashed hnes to the approximation 
from Eq. ([TT]). The curves are labeled by the value of e/3. 



Eq. (^ yields 
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with A = 3g — 2 — 2wa. As long as A > 0, the diffusion coefficients decrease with e. This 
happens when a < oq with oq = (3g — 2)/2w > a* > 0, i.e., when rj is either repulsive or 
slightly attractive relative to e. The decrease is again exponential: it behaves as exp{qe(3/2) 
for a < a* (the termp^ prevails) and as exp{Xef3/2) for a* < a < ao (the term exp{—wriP)p\^^ 
prevails). We illustrate this behavior in Fig. [7] for a = a* when both terms are important. 

However, if a > oq, then A < 0, and exp{—wr](3)pl,j^ always prevails over the term pg. In 
this case the diffusion coefficients grow exponentially fast with e as exp(A£/3/2). Therefore, 
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FIG. 8: The coverage dependences of the diffusion coefficients for the three lattices when rj = ae 
with a > qq. The soHd lines correspond to Eq. ([5]) and the dashed lines to the approximation from 
Eq. (jlip . The curves are labeled by the value of e/3. For the square and hexagonal lattices the 
approximation from Eq. (|11|) incorrectly increases with 9 near 9 = 0. 

an attractive e actually boosts diffusion as soon as rj is sufficiently attractive with respect to 
e — more than its oq multiple (see Fig. |8]). For the triangular lattice the threshold value ao 
is as large as 4, but for the square and hexagonal lattices it is comparable to 1, being equal 
to 5/4 and 7/8, respectively. 

D. Comparison with other results 

As mentioned in Section [TTl model ([T]) was previously used to study surface diffusion on a 
triangular and square lattices, using Monte Carlo simulations and real-space renormalization 
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FIG. 9: The chemical diffusion coefficient with the interaction in the activated state neglected 
{rj = 0) for (a) a triangular lattice and temperatures T = —Ke/kB with K = 1.35,1.14,1.0,0.94 



from top to bottom (adapted from Ref. 



15l): (b) a square lattice and 



with K = 0.63 (upper) and K = 0.57 (lower) (adapted from Ref. 



temperatures T = —Ke/kB 



14l). Squares represent Monte 



Carlo simulations, dashed lines a renormalization group method, and solid lines Eq. ([5]). 

group methods.—""— However, temperatures above the critical point were considered there, 
so a meaningful comparison with our low-temperatures results cannot be carried out. To 
illustrate this, in Fig. M we show the previously obtained results as well as our results for 
the same temperatures in the case when rj is neglected, boldly applying Eq. (|5]) in the 
over-critical region. Clearly, the agreement is only a qualitative and quickly worsens as the 
temperature grows. 

IV. CONCLUSIONS 

We investigated the interplay between the opposing effects of the interactions in the 
adsorbed and activated states on surface diffusion at subcritical temperatures. To this 
end, the coverage dependences of the chemical and jump diffusion coefficients D^ and Dj 
were considered and their sensitivity to the interactions was studied. We stressed that at 
low temperatures a phase transition is necessary whenever a coverage dependence of any 
quantity is to be explored; far from transitions a single phase is stable in the system and 
the coverage is almost constant. 

We considered a simple model in a finite portion of a triangular, square, and hexagonal 
lattice with periodic boundary conditions. An adsorbed particle was assumed to interact 
only with its q nearest neighbors via an attractive energy e, while an activated particle 
interacted with its w nearest adparticles via a different (attractive or repulsive) energy r]. 
The model is known to exhibit a first-order phase transition at low temperatures between 
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the fully vacant and fully occupied phase. Uncorrelated nn jumps were assumed. Thus, the 
activated-state interaction t] was associated with segments S oi w + 2 lattice sites, two of 
which represented a nn pair of vacant sites between which a jump was performed. 

For this model we were able to explicitly evaluate the dependence of Df. and Dj on the 
two interactions. It turned out that only two contributions were dominant, one for each 
phase: in the vacant phase it corresponded to the segment S in which all w + 2 sites were 
vacant; in the occupied phase it corresponded to S with the maximal number w of occupied 
sites. Three specific cases were discussed. 

First, if Tj was neglected, then, as expected, e decelerated diffusion. We showed that the 
rate of deceleration was exponential, as exp{qe/2kBT). 

Second, if e was fixed and r] varied, then an attractive rj accelerated diffusion only when it 
was more attractive than {q — l)/'w x e; again, the rate of acceleration was exponential, this 
time as exp{—wri/kBT). On the other hand, a less attractive or repulsive t] had practically 
no effect on diffusion. 

The most intriguing case occurred when the two interactions were considered to be propor- 
tional, 7] = ae. Then surface diffusion was exponentially accelerated (decelerated) whenever 
a > ao (a < ag), where the threshold value Oq = (3g — 2)/2w. Thus, an rj comparable to e 
is enough to boost diffusion on the square and hexagonal lattices (oq = 5/4 and oq = 7/8), 
while on the triangular lattice rj must be more than four times as attractive as e. 

The key ingredient underlying our investigation was Eq. ([5]). It allows one to study the 
case of any relation, r] = f{e), between the two interactions; the linear relation rj = ae 
was considered due to its simplicity. Moreover, since the equation is of rather general 
nature,"iS the approach presented here for model ([1]) can be applied also to other models on 
homogeneous or even heterogeneous lattices. This only requires to calculate the single-phase 
correlation factors and coverages for the studied model and at a given phase transition, which 
is simple to carry out heuristically, as was shown in Section IIIII A rigorous evaluation can 
be done via low-temperature cluster expansions.— 
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